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(54) Color printing method and apparatus using gamut mapping in munsell space. 

(57) A method and apparatus for constructing a 
color printer table. Printable colors in CIELAB 
space are transformed to Munsell space and 
color mapping functions are derived which map 
Munsell space coordinates into primary color 
coordinates CMY of the printer. A Munsell print- 
er gamut is determined by selecting points in 
Munsell space where the color mapping func- 
tions result within the printable primary color 
range. All discrete points of an extended gamut 
in CIELAB space are mapped into the Munsell 
space. The extended gamut comprises colors 
inside the printer gamut as well as colors out- 
side the printer gamut, such as colors typically 
found in a color monitor. Primary color values 
are inserted into the color printer table by tak- 
ing the mapped CIELAB points and applying the 
color mapping functions to the mapped CIELAB 
points that lie inside the Munsell printer gamut, 
and applying gamut mapping to the mapped 
CIELAB points in Munsell space that lie outside 
the Munsell printer gamut. Since the gamut 
mapping takes place in Munsell space, the hue 
is exactly preserved by simply preserving the 
hue angle. 
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The present invention pertains to a method and 
apparatus for building and using look-up tables which 
determine the colors that a color printer prints in re- 
sponse to requests to print specific colors. The spe- 
cific colors requested to be printed may include colors 5 
that are not printable by the printer. For those colors 
in particular, gamut mapping is performed in Munsell 
color space where hue planes are straight, i.e. not 
curved, so as to compensate for the Abney effect and 
to preserve the perceived hue. 10 

Recently, as the availability of color monitors and 
color printers has increased, it is more and more com- 
monplace for a computer user to view a full color im- 
age on a color monitor and then to request a full color 
printout of that image on a color printer. 15 

However, color printers and color monitors form 
color images differently. Specifically, a color monitor 
is a light emitting device; colors are formed on color 
monitors by adding light from three color primaries, 
generally, red, green and blue. Printed images, on the 20 
other hand, simply reflect ambient light; colors are 
perceived by the way ambient light is affected by 
three subtractive primaries, generally cyan, magenta ' 
and yellow (and sometimes black). 

These processes are fundamentally different. As 25 
a result, the range of colors displayable on the monitor 
is different from the range of colors printable by a 
printer. Figure 1 is the CIE 1931 chromaticity diagram 
showing the range (or "gamut") of colors displayable 
by a monitor (area "A") and the range (or "gamut") of 30 
colors printable by a printer (area "B"). As seen in Fig- 
ure 1, the range of colors displayable on a monitor is 
generally greater than the range of colors printable by 
a printer. This is because a monitor is a light emitting 
device and is able to display colors with greater sat- 35 
u ration. There are, however, some low saturation 
areas such as at area 10 where a printed image, 
which is light-subtractive, has greater color range 
than a monitor. 

Because of the difference between the ranges of 40 
printable and displayable colors, it has not heretofore 
been possible to print color images which are per- 
ceived as faithful reproductions of displayed color im- 
ages. Specifically, it is simply not possible to print a 
color in areas like out-of-gamut area 11 which are out- 45 
side the range "B" of printable colors. Accordingly, 
even though those colors may be seen on color mon- 
itors, they cannot be printed on a color printer. 

In U.S. Patent 4,941 ,038, out-of-gamut colors are 
adjusted to printable colors within the printergamut by 50 
selecting the printable color which has the shortest 
vector distance from the unprintable color and which 
preserves the hue angle of the unprintable color. Ex- 
periments and observations on human color percep- 
tion, however, have shown that a line representing a 55 
constant hue or color drawn from pure white out to the 
fully saturated color is not a straight line but rather is 
a curved line. The curvature of those lines (the so-cal- 



led "Abney Effect") is shown in the chromaticity dia- 
gram of Figure 2. For regions where the curvature of 
constant hue lines is relatively low, like bluish-green 
region 18, preserving hue angle only slightly changes 
the perceived hue. Thus, preserving hue angle for the 
bluish-green out-of-gamut color 18a by extending the 
hue back to point 18b on the printer gamut border 10 
still prints a bluish-green color. But for regions where 
the curvature of constant color lines is relatively high, 
like purplish-blue region 19, preserving hue angle 
greatly affects perceived hue. Thus, preserving hue 
angle for the blue out-of-gamut color 19a by extend- 
ing hue back to point 1 9b on the printer gamut border 
10 prints a color with a decidedly purplish-blue hue. 

Since each out-of-gamut color was adjusted inde- 
pendently, out-of-gamut colors are printed with poor 
color smoothness, where small changes in command- 
ed color can result in large changes in printed color. 
In particular circumstances, poor color smoothness 
manifests itself as non-monotonic changes in lumi- 
nance whereby the lightnesses of out-of-gamut colors 
does not increase smoothly and monotonically from 
dark to light but rather dips occasionally from light to 
dark. This results in a situation where colors which 
should merge. smoothly and monotonically from dark 
to light in fact show undesirable bands of darkness. 

In U.S. Patent 5,299,291 , which is assigned to the 
assignee of the present invention, to compensate for 
the Abney effect the hue angles in the printer table 
and border table are warped, and out-of-gamut colors 
are adjusted to printable colors within the printer ga- 
mut by building printer tables which produce colors 
which vary smoothly in the out-of-gamut regions and 
exhibited monotonic increases in lightnesses. Transi- 
tion colors are gamut mapped at a constant angle to 
colors of the same hue on the border of the printerga- 
mut. However, since the gamut mapping is performed 
in CIELAB space, it is not easy to see which colors are 
of the same hue because the hue planes are slightly 
curved (the Abney effect). This can be seen in Figure 
3, where the gradation of colors of the same hue from 
white to a fully saturated color are shown projected 
into the a* b* plane at a constant L*. The curved hue 
planes in CIELAB space are 5R (red), 5Y (yellow), 5G 
(green), 5B (blue) and 5P (purple), which are shown 
as curved lines. The warping computations in CIELAB 
space are inexact and time consuming. The contents 
of U.S. Patent 5,299,291 are incorporated herein by 
reference. 

A first aspect of the present invention addresses 
the foregoing difficulties and provides printer tables 
whose perceived hue is preserved and whose colors 
vary smoothly in the out-of-gamut regions by per- . 
forming gamut mapping in Munsell space where the 
hue planes are straight. Munsell space is shown in 
Figure 4 where hue planes 5R, 5Y, 5G, 5B and 5P are 
shown projected into the U W plane at a constant val- 
ue V as straight lines. Munsell coordinates are nor- 
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mally cylindrical coordinates H, V, C t but are shown in 
cartesian coordinates for a clearer understanding. 

The printer tables may be dervied once by a man- 
ufacturer, to be used in a printer driver which is sold 
to end users as part of a printer. 

In another aspect, the invention is a method and 
apparatus for building a printer table and a border ta- 
ble which give CMY values that are printed in re- 
sponse to a command to print a color. The CMY val- 
ues inserted into the respective tables are determined 
in Munsell space. Since the gamut mapping takes 
place in Munsell space, the hue is exactly preserved 
by simply preserving the hue angle, and out-of-gamut 
colors vary smoothly and exhibit monotonic increase 
in lightnesses for increasingly lighter out-of-gamut 
colors. 

This is accomplished by transforming measure- 
ments of printable colors to Munsell space. Color 
mapping functions are derived which map Munsell 
space coordinates into CMY coordinates. A Munsell 
printer gamut is determined by selecting points in 
Munsell space where the color mapping functions re- 
sult within the printable CMY color range. All discrete 
points of an extended gamut in CI E LAB space are 
mapped into the Munsell space. The extended gamut 
comprises colors within the printer gamut as well as 
colors outside the printer gamut such as colors typi- 
cally found in a color monitor. CMY values are insert- 
ed into the color printer table by taking the mapped 
CIELAB points and applying the color mapping func- 
tions to the mapped CIELAB points that lie inside the 
Munsell printer gamut, and applying gamut mapping 
to the mapped CIELAB points in Munsell space that 
lie outside the Munsell printer gamut. Since the gamut 
mapping takes place in Munsell space, the hue is ex- 
actly preserved by simply preserving the hue angle. 

If desired, a border table may be provided for col- 
ors outside the extended gamut by mapping border 
points of the extended gamut, at all discrete angles 
and lightnesses, into the Munsell space and applying 
gamut mapping to the mapped CIELAB points in Mun- 
sell space. 

In another aspect, color printing is conducted by 
reference to the printer table and the border table, 
which give CMY values in response to a command to 
print a color. The printer table gives CMY values for 
colors within a printer table domain which consists of 
rectangular areas surrounding the extended gamut. 
The border table gives CMY values for colors outside 
the printer table domain. 

This brief summary has been provided so that the 
nature of the invention may be understood quickly. A 
more complete understanding of the invention can be 
obtained by. reference to the following detailed de- 
scription of the preferred embodiment thereof in con- 
nection with the drawings. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a chromaticity diagram showing how 
the gamut of colors printable on a printer is relat- 
5 ed to the gamut of colors displayable on a moni- 

tor; 

Figure 2 is a chromaticity diagram showing the 
curved constant hue lines that characterize the 
Abney effect; 

10 Figure 3 is a graph illustrating the gradation of 

colors of the same hue in CIELAB space; 
Figure 4 is a graph illustrating the gradation of 
colors of the same hue in Munsell space; 
Figure 5 is a block diagram of a printing appara- 

15 tus according to the invention; 

Figure 6 is a flow diagram showing how a printer 
driver in the Figure 6 apparatus selects CMYK 
values for a color printer; 

Figures 7(a) - 7(d) when taken together as shown 
20 in Figure 7, comprise a flow diagram for describ- 

ing how the printer table and the border table are 
constructed; 

Figure 8 shows how unconnected regions are re- 
moved from the Munsell printer table; 
25 Figure 9 shows radially concave regions in the 

Munsell printer gamut; 

Figure 10 illustrates how to remove concavities 
by making the Munsell printer table radially con- 
vex; 

30 Figure 11 illustrates how the border points of the 

extended gamut in CIELAB space are chosen to 
determine the border table; 
Figure 12 shows the gamut mapping in the Mun- 
sell space; 

35 Figure 1 3 illustrates a typical division of CIELAB 

space into a printer table; 

Figure 14 illustrates the arrangement of a border 

table in CIELAB space; and 

Figure 15 shows the relation between the printer 
40 gamut, the extended gamut, the rectangular area 

of the printer table and the wheel-shaped border 

table in CIELAB space for an arbitrary luminance 

value L*. 

45 DETAILED DESCRIPTION OF THE PREFERRED 
EMBODIMENT 

Figure 5 is a block diagram of a printing appara- 
tus according to the invention. 

50 As shown in Figure 5, the printing apparatus in- 

cludes a host CPU 20, a color monitor 30 and a color 
printer 40. Host CPU 20 includes a processor 21 such 
as an 80586 microprocessor, a random access mem- 
ory ("RAM") 22 which provides working storage area 

55 to processor 21, a read only memory ("ROM") 24 
which provides static storage for processor 21, mon- 
itor driver 25 and a printer driver 26. Host CPU 20 is 
accessed by an operator via keyboard 27 which is 
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connected through an interface 29 to processor 21. 
Using the keyboard, an operator can cause processor 
21 to execute stored program instructions which 
cause color images to be displayed on monitor 30 and 
which cause corresponding color images to be print- 5 
ed on color printer 40. 

Other peripheral devices, such as disk drives, 
tape drives, color video interfaces, color scanner in- 
terfaces, etc., may be provided for host CPU 20 but 
those other devices are not shown in the interest of 10 
simplicity. In cooperation with the stored program in- 
structions executed by processor 21, such devices 
permit, for example, a color image to be scanned into 
RAM 22 and displayed on monitor 30, the colors in the 
image to be manipulated, and the resulting image to 15 
be printed on printer 40. 

In accordance with stored program instructions, 
processor 21 derives a color image for display on 
monitor 30. Processor 21 provides the color image to 
monitor driver 25 which in turn derives RGB values for 20 
each pixel in monitor 30. The RGB values are provid- 
ed via interface 31 to the monitor 30 where those val- 
ues are displayed. 

Upon request, processor 21 also feeds a color im- 
age to printer driver 26 for printing by color printer 40. 25 
Printer driver 26 derives CMY values for each pixel of 
the color image based on the color values provided 
from processor 21. The CMY values are determined 
in accordance with either the printer table 26a or the 
border table 26b. The printer table 26a is a table which 30 
provides CMY values for all colors that are printable 
by printer 40. The border table 26b is a table which 
provides suitable CMY values for colors that are not 
printable by printer 40. The printer table 26a may also 
include CMY values for some unprintable colors so as 35 
to smooth the transition from printable to unprintable 
colors. That is, printer table 26a has a domain which 
includes an extended gamut which includes the prin- 
ter gamut and a typical color monitor gamut. In addi- 
tion, a black (hereinafter "K") value may also be de- 40 
rived. The CMYK values are fed via interface 41 to 
printer 40 where they are stored in bit map memory 
42 within printer 40. The bit map memory 42 may 
store a full bit map image of the printed image or it 
may store only a band or partial bit map image. When 45 
sufficient color data is stored in bit map memory 42, 
a color printer head 44 reciprocates across a platen 
adjacent to a sheet of paper. In the present embodi- 
ment, print head 44 includes 32 ink jet nozzles which 
may be arranged in a four column by eight row pat- 50 
tern. The nozzles in the first column all eject droplets 
of cyan ink; the nozzles in the second column all eject 
droplets of magenta ink; the nozzles in the third col- 
umn all eject droplets of yellow ink; and the nozzles 
in the fourth column all eject droplets of black ink. The 55 
nozzles are controlled independently in accordance 
with the color data in bit map memory 42 such that in 
one reciprocation of print head 44 across the platen, 
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eight rows of pixels are printed. 

Figure 6 is a flow diagram showing how printer 
driver 26 selects CMYK values from the color data 
provided by processor 21 . In step S601 , printer driver 
26 receives RGB values for a location (x,y) in bit map 
memory 42. In step S602, printer driver 26 drives de- 
vice independent color coordinates from the RGB val- 
ue. Preferably, the device independent coordinates 
are CIELAB coordinates. This is because the coordin- 
ates in CIELAB space are perceptually uniform such 
that equal-sized intervals anywhere in CIELAB space 
correspond to equal-sized changes in perceived col- 
or. Furthermore, CI ELAB space is one of the standard 
color spaces utilized in the printing industry. 

In step S603, the luminance coordinate is com- 
pressed at the extremes of the L* axis in CIELAB 
space. Compression step S603 may be performed 
explicitly by mathematical manipulation of the L* val- 
ue from step S602 or implicitly by storing modified 
CMY values in the printer table and the border table. 
If performed implicitly, which is preferable in some in- 
stances, then both the printer table and the border ta- 
ble store pre-compressed values. More particularly, 
the printer table and the border table can be arranged 
so that the values stored at, for example, luminance 
L* = 99 actually correspond to a luminance of L* = 94. 
Likewise, values stored at, say, luminance L* = 7 ac- 
tually correspond to a luminance of L* = 26. Values in 
the center of the luminance range, from for example 
L* = 38 through 90, remain unmodified. This arrange- 
ment results in luminance compression without the 
need for explicit compression. 

While compression step S603 is optional, it is 
nevertheless preferable to perform since it ensures 
that colors at extreme values of luminance are printed 
with perceptible changes in luminance. More particu- 
larly, because monitor 30 displays colors with light 
emitting elements, it can display colors with higher 
values of luminance then those achievable by printer 
40, whose highest value of luminance is limited by the 
whiteness of the paper upon which the color image is 
formed. Further still, since monitor 30 can completely 
turn off its light emitting elements, it can display col- 
ors with lower values of luminance than those print- 
able by printer 40, since even black ink reflects some 
ambient light. Accordingly, to ensure that some color 
is printed, even at the highest and lowest luminance 
values, it is preferable to compress the luminance val- 
ues determined in step S602 into a range that is print- 
able by printer 40. 

in step S604, the L* f a* and b* coordinates de- 
rived in steps S602 and S603 are inspected to deter- 
mine whether they fall within the range covered by 
printer table 26a. If the L*. a*, b* coordinates are with- 
in the range covered by printer table 26a, then flow 
advances to step S605 which looks up the corre- 
sponding CMY values in printer table 26a at location 
L*, a*, b* (actually, the nearest location since only dis- 
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crete values of L*, a* and b* are stored). On the other 
hand, if the L*, a*, b* coordinates are not within the 
printer table 26a, then flow advances to step S606 in 
which the hue angle 9 is derived from the a* and b* 
values according to the following formula: 5 
9 = arctan (b * /a * ) 

Flow then advances to step S607 which looks up 
corresponding CMY values in border table 26b at the 
nearest location which corresponds to the luminance 
L* and the hue angle derived in step S606. 10 

In either event, flow then advances to step S608 
in which the CMY values are stored in bit map mem- 
ory 42 at location (x,y). If desired, the CMY values 
may be modified before storage, for example, by in- 
terpolation, so as to accommodate the difference be- 15 
tween the actual L*, a*, b* values stored in the tables 
and the desired values calculated above. 

In step S609, printer driver 26 determines wheth- 
er the bit map memory has been completed. If the bit 
map memory has not been completed, then flow re- 20 
turns to step S601 in which the next RGB value is re- 
ceived for the next location (x,y) in bit map memory. 
On the other hand, if the bit map memory has been 
completed, or if a sufficient area of the bit map mem- 
ory has been completed (such as an eight row long 25 
band corresponding to the eight rows of ink jet noz- 
zles in head 44), then flow advances to step S610 
where gamma correction is performed. Gamma cor- 
rection adjusts the CMY values in the bit map memory 
so as to achieve a uniform distribution of luminance. 30 
In step S611, undercolor removal in the present em- 
bodiment may be performed by the simple expedient 
of selecting the minimum value of CMY and assigning 
that value to the black value. Then, each of the CMY 
values is adjusted by subtracting the black value from 35 
it. 

The order of steps S610 and S611 is not critical 
and those steps may be switched, for example to ac- 
commodate a particular color printing technique such 
as continuous tone, dither techniques or error diffu- 40 
sion. 

In step S612, color printing is initiated using the 
resulting CMYK values. 

Figures 7(a) - 7(d) when taken together as shown 
in Figure 7, comprise a flow diagram for describing 45 
how the printer table 26a and the border table 26b is 
constructed. The flow procedures shown in Figures 
7(a) - 7(d) need only be performed once to derive prin- 
ter tables and border tables for each printer. The prin- 
ter tables are derived once by a manufacturer, to be so 
used in a printer driver which is sold to end users as 
part of a printer. The printer tables may be in the form 
of software and/or hardware. 

The printer and border tables are CIELAB printer 
tables and CIELAB border tables because the RGB 55 
values received from an input device such as a scan- 
ner are generally transformed into CIELAB coordin- 
ates. However, the CMY values which are inserted 



into the CIELAB printer table and the CIELAB border 
table are computed in the Munsell space where hue 
planes are straight, to exactly preserve the hue by 
merely preserving the hue angle. 

The steps shown in the flow diagrams of Figure 
8 may be stored program instructions which are op- 
erable in a CPU such as the CPU 20 shown in Figure 
5, or the steps may be executed entirely in a special 
purpose hardware apparatus. 

In step S701, relative spectral reflectance meas- 
urements are performed on the range of colors print- 
able by printer 40. Preferably, this is achieved by 
printing a very large subset, or a complete set, of all 
colors printable by printer 40. For example, in the prin- 
ter used in the present embodiment, each of the CMY 
and K values may be printed in 65 gradations ranging 
numerically from 0 to 64. A subset of about one quar- 
ter of those values, for each color, are printed. Thus, 
for example, 17 C values are printed, namely numer- 
ical values 0, 4, 8, 12, ... 64, and 17 M values are print- 
ed, and 17 Y values are printed. All possible combin- 
ations of those 17 CMY values are printed, yielding 17 
x 17 x 17 = 4,913 color patches. 

In addition to the foregoing colors, all possible 
gray values, in this case 48 additional gray values 
over the 17 already printed, are also printed. 

With the foregoing sampling of the printer gamut, 
it will be seen that pure gray colors are printed togeth- 
er with hued colors. Whatever method of sampling is 
used, this property of pure gray printing should be 
preserved since proper gray reproduction is a desir- 
able property in color reproduction. 

Relative spectral reflectance measurements are 
performed on each of the 4,91 3 color patches and 48 
additional gray patches. 

Step S702 transforms the measurements of 
printable colors into points in the CIELAB space with 
respect to CIE standard illuminant C in L*, a* and b* 
coordinates which are transformed into H, V, and C 
cylindrical coordinates in Munsell space in step S703. 

In step S704, the cylindrical H, V, C coordinates 
in Munsell space are transformed to U, V, W cartesian 
coordinates, where the U, W plane is orthogonal to 
the V axis, with U pointing in the direction of the Mun- 
sell hue 1 0RP and W pointing in the direction of Mun- 
sell hue 5Y. 

Step S705 derives color mapping functions which 
map the Munsell coordinates into primary color coor- 
dinates, such as CMY coordinates. In the present em- 
bodiment, a cubic least square fit from Munsell space 
into CMY space was chosen. That is, using well- 
known least square fitting techniques, coefficients Cq 
through c 19 , m 0 through m 19 , and y 0 through y 19 were 
derived to give the best fit, in the least squared sense, 
to the printable colors measured in step S701 and 
transformed in step S702: 
C = Cq + c,V + C2U + c 3 W + C4V 2 + C5U 2 + 

CsW 2 
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+ c 7 VU + c 8 VW + CgUW + c 10 V 3 + c^U 3 + 

c 12 W 3 + c 13 V2U + c 14 VU 2 + c 15 V2W + 
c 18 VW2 + c 17 U2W + c 18 UW2+ c 19 VUW (1) 
M = m 0 + + m 2 U + m 3 W + m 4 V 2 + m 5 U 2 
+ m 6 W 2 

+ m 7 VU + m 8 VW + m 9 UW + m 10 V 3 + m^U 3 + 

m^W 3 + m^U + m 14 VU 2 + m^V^W + 
m 16 VW 2 + m 17 U 2 W + m 18 UW 2 + m 19 VUW (2) 
Y = y 0 + YiV + y 2 U + y 3 W + y 4 V 2 + y 5 U 2 + 

y 6 W2 

+ y 7 VU + y 8 VW + y 9 UW + y 10 V 3 + y^U 3 + 
yi 2 W 3 + y^U + y 14 VU 2 + y 15 V 2 W + 

yieVW 2 + y 17 U 2 W + y^UW 2 + y 19 VUW (3). 
In step S705, any mathematical function which 
fits the measurements taken in step S801 from the 
device independent coordinate space to CMY coor- 
dinate space may be used. Preferably, however, the 
mapping function includes smoothing so as to elim- 
inate measurement irregularities that may have been 
encountered in step S701. 

It may, in addition, be preferable to weight some 
of the points measured in step S701 prior to deriving 
mapping in step S705. For example, proper skin tone 
color reproduction is an important property of color 
printers. Accordingly, it may be desirable, in some cir- 
cumstances, to weight colors in the area of skin tone 
colors more heavily than other colors. 

In step S706, the device independent space, 
namely Mgnsell space, is divided into equally sized 
discrete intervals, one of the intervals including the V 
axis such as by being centered at the V axis. 

The determining of the Munsell printer gamut is 
started in step S707 by determining the points where 
the color mapping functions result within the printable 
primary color range, i.e., the CMY values. The size of 
the intervals in the Munsell printer gamut should be 
made as small as possible. For example, it has been 
found that fine luminance gradations are more impor- 
tant than fine hue and saturation gradations. It has 
been determined that dividing the luminance axis into 
intervals of AV = 0.1 (value V ranges from 0 through 
10) provides adequate value V gradation. On the 
other hand, such fine gradations are not ordinarily 
needed in hue and saturation, and therefore AU = AW 
= 0.4 has been found to yield adequate hue and sat- 
uration gradations (U and W range from about -20 
through +20 near the center (V = 5) of the value V 
axis). 

In addition to the foregoing considerations, it 
should also be observed that the Munsell printer ga- 
mut is not the same for each value V. Specifically, the 
gamut is relatively smaller at value V extremes and 
relatively larger at the center of the value axis. 

In steps S708 and S709, the Munsell printer ga- 
mut is adjusted. In Figure 8 it is shown that the color 
matching function used may give rise to false regions 
in the Munsell printer gamut, such as region 45, which 
is not within the true Munsell-printer gamut 46. These 



artifacts are removed in step S708 by removing all re- 
gions that are not connected to the region around the 
V axis. 

In Figure 9 the Munsell printer gamut shown is not 
5 radially convex because each and every radial line 
from the V axis does not intersect border 47 at one 
and only one point. In particular, radial line 48 inter- 
sects border47 at three points 49a, 49b and 49c. The 
region between points 49a and 49b is a radial concav- 
10 ity and can cause the generation of inappropriate 
CMY values. Accordingly, in step S709, the values in 
the Munsell printer gamut are adjusted to make them 
radially convex. 

Figure 10 illustrates this process. Figure 10 
15 shows the rectangular grid in the U and W axis for an 
arbitrary luminance value V. Cells 51 through 55 all in- 
clude printable values within the Munsell printer ga- 
mut. However, cell 57 is a radially concave cell be- 
cause the color mapping functions result outside the 
20 primary color range in values C = 9, M = 11 and Y= - 
1, and the radial line at angle H crosses two cells in 
the Munsell printer gamut (cells 53 and 55). Accord- 
ingly, a CMY value is assigned to cell 57 to make the 
table radially convex. The CMY value is selected by 
25 simply truncating the result to C = 9, M = 11 and Y = 
0, which places cell 57 inside the primary color range. 

The printer table 26a and border table 26b pro- 
vide CMY values in a device independent space, 
namely the CIELAB space. Prior to later insertion of 
30 CMY values into these tables, the CIELAB space 
must be divided into equally sized discrete intervals, 
one of the intervals including the L* axis such as by 
being centered at the L* axis. In step S710 the CIE- 
LAB space is so divided. 
35 In step S711, the extended gamut in the CIELAB 

space is determined by first determining the printer 
gamut in the CIELAB space and then selecting the ex- 
tended gamut to be slightly larger than the printer ga- 
mut and a typical color monitor gamut. The printer ga- 
40 mut in the CIELAB space is determined by applying 
the color mapping functions which map the CIELAB 
coordinates into primary color coordinates, i.e. CMY 
coordinates, and making the CIELAB printer gamut 
connected to the L* axis and radially convex. 
45 In step S712 all discrete points of the extended 

gamut in the CIELAB space are mapped into the Mun- 
sell Space. In step S713 the border points of the ex- 
tended gamut in the CIELAB space are mapped at all 
discrete angles and lightnesses into the Munsell 
so space. Figure 11 shows the border points of the ex- 
tended gamut 80 in the CIELAB space. In areas like 
the area 81 at the border of the extended gamut 80 
between points 82 and 83, the saturation changes 
too much. To avoid this large saturation change when 
55 going from one discrete angle to the next, sufficiently 
fine angle changes have to be chosen. Increments of 
1° have been found to be sufficiently fine. 

In step S714 the size of the extended gamut in the 



6 



11 



EP 0 682 440 A2 



12 



CIELAB space is reduced, if necessary, so that the 
mapped points in the Munsell space are clearly de- 
fined Munsell colors, which means that they lie within 
the area of the Munsell color patches, or within the 
area of extrapolated Munsell colors. The extrapola- 5 
tion of Munsell colors is described in the "Final Report 
On The O.S.A. Subcommittee On The Spacing Of 
Munsell Colors" by Newhall, S. et al., J. Opt. Soc. Am., 
33,7(1943). 

At step S715 primary color values, i.e. CMY val- 10 
ues are inserted into the CIELAB printer table around 
the L* axis by taking the mapped CI ELAB printer table 
points and applying the color mapping functions to 
the mapped CIELAB points that lie inside the Munsell 
printer gamut, and applying gamut mapping to the 15 
mapped CIELAB points in Munsell space that lie out- 
side the Munsell printer gamut. The gamut mapping 
maps a color point in the Munsell space which lies 
outside the Munsell printer gamut at a constant angle 
to a color of the same hue at the border of the Munsell 20 
printer gamut. Since gamut mapping takes place in 
the Munsell space, the hue is exactly preserved by 
simply preserving the hue angle. 

How gamut mapping is performed may be seen 
by referring to Figure 12. In Figure 12 a hue plane in 25 
the Munsell space is shown, and 60 is the border of 
the printer gamut. Each mapped point in the Munsell 
space that lies outside the Munsell printer gamut 60, 
such as the point 51 is gamut mapped at a constant 
angle «: to a color of the same hue at a point on the 30 
border 60, such as point 52. The gamut mapping as- 
signs to point 51 the CMY color of point 52, which is 
determined by applying the color mapping functions 
to point 52. Likewise, a mapped point 53 outside the 
Munsell printer gamut 60 is gamut mapped at the con- 35 
stant angle oc to a color of the same hue at point 54 
on the border 60. The gamut mapping assigns to point 
53 the CMY color of point 54, which is determined by 
applying the color mapping functions to point 54. The 
angle oc lies in the range of approximately 1 0° to 30° 40 
with respect to the lightnesses plane at the points 51 
or 53. 

Within the displayed hue plane, point 50 denotes 
the point of maximum saturation at the border 60 of 
the Munsell printer gamut. For points above the point 45 
50, such as the point 51, the angle oc extends down- 
wardly; and for points below the point 50, such as the 
point 53, the angle oc extends upwardly. All mapped, 
points in the wedge 70, subtended by the angle 2oc 
from the point 50 are gamut mapped to the color at 50 
point 50. The gamut mapping assigns to all points in 
the wedge 70 the CMY color of point 50, which is de- 
termined by applying the color mapping functions to 
point 50. Since the gamut mapping preserves the hue 
angle in Munsell space, the hue is exactly preserved. 55 

At Step S716, CMY values are inserted into the 
printer table at the L* axis by determining the printer 
grays. 



At step S717, CMY values are inserted into the 
border table by taking the mapped CIELAB border 
points of the extended gamut and applying the gamut 
mapping as shown in Figure 1 1 . The border table lists 
one primary CMY color at all discrete angles and light- 
nesses for points outside the extended gamut in CIE- 
LAB space. These primary colors are determined by 
mapping the border points of the extended gamut in 
CIELAB space at all discrete angles and lightnesses 
into the Munsell space and applying the gamut map- 
ping in the Munsell space as described relative to Fig- 
ure 12. 

Figure 15 shows the correspondence of border 
tables and printer tables. Whereas printer table 26a is 
a rectangular grid for an arbitrary luminance value L* f 
border table 26b is arranged as a wheel centered at 
a* = b* = 0. Individual cells in the border table are ac- 
cessed by the angle 9 in the a* b* plane which also 
corresponds to hue. Experimentally, it has been 
found that 1 degree increments, resulting in 360 cells 
in each border table, provides adequate gradation of 
hue. However, this gradation may be adjusted. 

In step S718, the printer table is rectangularized. 
More specifically, until this step, CMY values have 
been inserted into the printer table only in areas within 
the extended gamut 61 which includes the printer ga- 
mut 60 (steps S710 and S715). In step S718, the re- 
maining cells of the printer table such as cells like cell 
68 in Figure 1 5 are filled out by calculating the hue an- 
gle for each blank cell that remains in the printer table 
and by inserting the border table color at that hue an- 
gle as illustrated at 69 in Figure 15. 

Whereas the printer table was arranged as a rec- 
tangular grid in the a* b* axes for each luminance val- 
ue, the border table is arranged in wheels with one 
wheel for each of the luminance values of the printer 
table. Thus, as shown in Figure 14, a wheel-shaped 
border table is provided for each of the luminance val- 
ues for which a printer table exists, which provides 
one border table in correspondence with each of the 
printer tables. The border tables contain plural cells 
which are accessed based on hue angle 9 calculated 
as a function of the a* and b* coordinates, as follows: 
9 = arctan (b */a * ) 

Figure 13 illustrates a typical division of CIELAB 
space into a printer table, although all luminance and 
hue gradations have not been shown to simplify the 
presentation. At relatively low luminance values such 
as L* = 10, a relatively small rectangular grid in the a* 
and b* axes is adequate to map the printer gamut. 
Similarly, at relatively high luminance values, such as 
L* = 90, a relatively small rectangulargrid in the a* and 
b* axes is also adequate to store the printer gamut. 
However, at intermediate luminance values, such as 
that at L* = 50, a relatively larger rectangular grid in 
the a* and b* axes is required to map the printer ga- 
mut. 

As further shown in Figure 13, the rectangular 
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grid at each luminance level Includes the L* axis (In 
Figure 13 it Is centered on the L* axis). That is, there 
is a cell in the rectangular grid that corresponds ex- 
actly to a* = b* = 0. That central point, namely a* = b* 
= 0, corresponds to a pure gray color which, as men- 
tioned above, is desirably reproduced as a pure gray 
color for proper color reproduction. 



Claims 
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1. A method of constructing a printer table, compris- 
ing the steps of: 

transforming measurements of printable 
colors into Munsell space; 15 

deriving color mapping functions which 
map Munsell space coordinates into primary col- 
or coordinates; 

determining a Munsell printer gamut by se- 
lecting points in Munsell space where the color 20 
mapping functions result within the printable pri- 
mary color range; 

determining a CIELAB space extended 
gamut which includes a printer gamut and a typ- 
ical color monitor gamut; 25 

mapping all discrete points of the extend- 
ed gamut in CIELAB space into the Munsell 
space; and 

inserting primary color values into said 
printer table by gamut mapping the mapped CiE- 30 
LAB points in Munsell space that lie outside the 
Munsell printer gamut, at a constant angle to a 
color of the same hue at the printer gamut border 
in Munsell space. 

35 

2. The method of Claim 1, including the step of: 

inserting additional primary color values 
into said printer table by taking the mapped CIE- 
LAB points and applying the color mapping func- 
tions to the mapped CIELAB points that lie inside 40 
the Munsell printer gamut. 

3. The method of Claim 2, including the steps of: 

mapping border points of the extended ga- 
mut in CIELAB space into the Munsell space; and 45 

inserting primary color values into a bor- 
der table by applying the gamut mapping in Mun- 
sell space to the mapped CIELAB border points. 

4. A method of constructing a printer table, compris- so 
ing the steps of: 

transforming relative spectral reflectance 
measurements of printable colors into Munsell 
space; 

deriving color mapping functions which 55 
map Munsell space coordinates into primary col- 
or coordinates; 

determining a Munsell printer gamut by se- 

8 



lecting points in Munsell space where the color 
mapping functions result within the printable pri- 
mary color range; 

determining a CIELAB space extended 
gamut which includes the printergamut and a typ- 
ical color monitor gamut; 

mapping all discrete points of the extend- 
ed gamut in CIELAB space into the Munsell 
space; and 

inserting primary color values into said 
printer table by taking the mapped CIELAB points 
and applying the color mapping functions to the 
mapped CIELAB points that lie inside the Munsell 
printer gamut, and applying gamut mapping to the 
mapped CIELAB points in Munsell space that lie 
outside the Munsell printer gamut, 

wherein the gamut mapping maps a color 
point, in the Munsell space which lies outside the 
Munsell printer gamut, at a constant angle to a 
color of the same hue at the border of the Munsell 
printer gamut. 

5. The method of claim 4, including the steps of: 

mapping border points of the extended ga- 
mut in CIELAB space, at all discrete angles and 
lightnesses, into the Munsell space; and 

inserting primary color values into a bor- 
der table by applying gamut mapping in Munsell 
space to the mapped CIELAB border points. 

6. The method of Claim 4, including the step of: 

transforming the points in- CIELAB space 
into points in Munsell space for each printable 
color. 

7. The method of claim 6, including the step of: 

transforming the points in Munsell space 
in cylindrical coordiantes H, V, C into Cartesian 
coordinates U, V, W for each printable color. ■ 

8. The method of claim 7, including the step of: 

dividing the Munsell space into discrete in- 
tervals centered on the V axis. 

9. The method of claim 8, including the step of: 

adjusting the Munsell printergamut by re- 
moving unconnected regions from the Munsell 
printer gamut. 

10. The method of claim 9, including the step of: 

adjusting the Munsell printer gamut by 
making its radially convex. 

11. The method of claim 10, including the step of: 

dividing CIELAB space into discrete inter- 
vals centered on the L* axis. 

12. The method of claim 11, including the step of: 
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reducing the size of the extended gamut in 
the CIELAB space, if necessary, so that the map- 
ped points in Munsell space are clearly defined 
Munsell colors. 

13. The method of claim 12, including the step of: 

inserting primary color values into the 
printer table at the L* axis by determining printer 
grays. 

14. The method of claim 13, including the step of: 

inserting primary color values into a bor- 
der table by taking the mapped CIELAB border 
points and applying gamut mapping. 

15. The method of claim 14, including the step of: 

rectangularizing the printer table based on 
the border table. 

16. An apparatus for constructing a printer table, 
comprising: 

means for transforming measurements of 
printable colors into Munsell space; 

means for deriving color mapping func- 
tions which map Munsell space coordinates into 
primary color coordinates; 

means for determining a Munsell printer 
gamut by selecting points in Munsell space 
where the color mapping functions result within 
the printable primary color range; 

means for determining a CIELAB space 
extended gamut which includes a printer gamut 
and a typical color monitor gamut; 

means for mapping all discrete points of 
the extended gamut in CIELAB space into the 
Munsell space; and 

means for inserting primary color values 
into said printer table by gamut mapping the map- 
ped CIELAB points in Munsell space that lie out- 
side the Munsell printer gamut, at a constant an- 
gle to a color of the same hue at the printer gamut 
border in Munsell space. 

17. The apparatus of Claim 16, including: 

means for inserting additional primary col- 
or values into said printer table by taking the map- 
ped CIELAB points and applying the color map- 
ping functions to the mapped CIELAB points that 
lie inside the Munsell printer gamut. 

18. The apparatus of Claim 17, including: 

means for mapping border points of the 
extended gamut in CIELAB space into the Mun- 
sell space; and 

means for inserting primary color values 
into a border table by applying the gamut map- 
ping in Munsell space to the mapped CI ELAB bor- 
der points. 



19. Apparatus for constructing a printer table, com- 
prising: 

means for transforming relative spectral 
reflectance measurements of printable colors 
5 into Munsell space; 

means for deriving color mapping func- 
tions which map Munsell space coordinates into 
primary color coordinates; 

means for determining a Munsell printer 
10 gamut by selecting points in Munsell space 

where the color mapping functions result within 
the printable primary color range; 

determining a CIELAB space extended 
gamut which includes the printergamut and a typ- 
15 ical color monitor gamut; 

means for mapping all discrete points of 
the extended gamut in CIELAB space into the 
Munsell space; and 

means for inserting primary color values 
20 into said printer table by taking the mapped CIE- 

LAB points and applying the color mapping func- 
tions to the mapped CIELAB points that lie inside 
the Munsell printer gamut, and applying gamut 
mapping to the mapped CIELAB points in Mun- 
25 sell space that lie outside the Munsell printer ga- 

mut, 

wherein the gamut mapping maps a color 
point, in the Munsell space which lies outside the 
Munsell printer gamut, at a constant angle to a 
30 color of the same hue at the border of the Munsell 

printer gamut. 

20. The apparatus of claim 19, including: 

means for mapping border points of the 
35 extended gamut in CIELAB space, at all discrete 

angles and lightnesses, into the Munsell space; 
and 

means for inserting primary color values 
into a border table by applying gamut mapping in 
40 Munsell space to the mapped CIELAB border 

points. 

21. The apparatus of Claim 19, including: 

means for transforming the points in CIE- 
45 LAB space into points in Munsell space for each 

printable color. 

22. The apparatus of Claim 21, including: 

means for transforming the points in Mun- 
50 sell space in cylindrical coordinates H, V, C into 

Cartesian coordinates U, V, W for each printable 
color. 

23. The apparatus of claim 22, further comprising: 
55 means for dividing the Munsell space into 

discrete intervals centered on the V axis. 

24. The apparatus of claim 23, further comprising: 
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means for adjusting theMunsell printer ga- 
mut by removing unconnected regions from the 
Munsell printer gamut 

25. The apparatus of claim 24, further comprising: 5 

means for adjusting the Munsell printer by 
making it radially convex. 

26. The apparatus of claim 25, further comprising: 

means for dividing CI E LAB space into dis- 10 
crete intervals centered on the L* axis. 

27. The apparatus of claim 26, further comprising: 

means for reducing the size of the extend- 
ed gamut in the CIELAB space, if necessary, so 1 5 
that the mapped points in Munsell space are 
clearly defined Munsell colors. 

28. The apparatus of claim 27, further comprising: 

means for inserting primary color values 20 
into the printer table at the L* axis by determining 
printer grays. 

29. The apparatus of claim 28, further comprising: 

means for inserting primary color values 25 
into a border table by taking the mapped CIELAB 
border points and applying gamut mapping. 

30. The apparatus of claim 29, further comprising: 

means for rectangularizing the printer ta- 30 
ble based on the border table. 

31. A method or apparatus according to any preced- 
ing claim wherein the constructed printer table is 
supplied in a package with a corresponding prin- 35 
ter apparatus. 

32. A colour printer package including a colour printer 
apparatus and means conveying printer table 
constructed in accordance with any preceding 40 
claim. 

33. A method or apparatus having the features of any 
combination of the preceding claims. 

45 
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